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DESCRIPTION 

Apparatus and method for the processing of sectional images 

The invention relates to a method and an apparatus for the processing of sectional images 
that are reconstructed from X-ray projections of an object, preferably from cone beam 
5 projections. 

Reconstructed 3D sectional images in C-arm based volume imaging often exhibit spatially 
slowly varying inhomogeneities caused by inconsistent projection data due to various 
reasons. First, scattered radiation constitutes one of the main problems. Especially for 
10 system geometries with large cone angle and therefore a large irradiated area, such as C-arm 
based volume imaging, scattered radiation produces a significant background added to the 
desired detected signal and thus significant cupping shaped inhomogeneities occur in the 
reconstructed volume. Further sources of inhomogeneities are beam hardening (cupping) 
and insufficient truncation correction (cupping and inverse cupping). 

15 

When a tight gray level window is applied to the viewed slices in the reconstructed 3D image 
in order to find low contrast details, the inhomogeneities cause parts of the region of interest to 
be either below or above the edges of the applied window. Low contrast details are therefore 
much more difficult to find in presence of inhomogeneities. 

20 

Anti-scatter grids reduce the amount of detected scattered radiation, but the reduction is 
not sufficient for typical system geometries for volume imaging, by additionally increasing 
the signal-to-noise ratio. Different algorithms for projection based scatter compensation 
have been proposed (cf. Maher K. P., Malone J. F., "Computerized scatter correction in 
25 diagnostic radiology", Contemporary Physics, vol. 38, no.2, pp. 131-148, 1997), but 
accurate quantitative scatter estimation is difficult to achieve. Truncation correction 
schemes like elliptical extension of projections (as described in R. M. Lewitt, Processing of 
incomplete measurement data in computed tomography, Med. Phys., vol. 6, no. 5, pp. 
412-417, 1979) can prevent severe shading artifacts at the boundaries of the FOV, but do 
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not completely correct for global shading. Beam hardening correction is possible, but 
requires a deep knowledge of acquisition parameters, that is not always available. 

Based on this situation it was an object of the present invention to provide simple means 
5 for improving sectional X-ray images that are affected by inhomogeneities. 

This object is achieved by an apparatus with the features of claim 1 and a method with the 
features of claim 10. Preferred embodiments are subject of the dependent claims. 

10 The apparatus according to the present invention is suited for the processing of sectional 

images of an object, wherein such sectional images are reconstructed from X-ray projections 
of the object that were taken from different directions. The apparatus may optionally be 
adapted to execute the reconstruction of the sectional images from projections or it may 
alternatively receive already reconstructed sectional images for further processing. The 

15 apparatus is adapted to execute the following steps with a sectional image tinder 
consideration: 

The determination of a baseline function that at least approximately describes 
spatially slowly varying artifacts of the image. As will be explained in more detail 
20 below, such kind of artifacts are quite typical for a class of sectional images 

reconstructed from X-ray projections. 

Calculating a corrected sectional image by compensating the original sectional 
image with the help of said baseline function, e.g. by subtracting the baseline 
25 function from the original image. 

The apparatus has the advantage to (approximately) correct certain slowly varying artifacts, 
wherein the correction particularly makes the inspection of a small gray-value window 
possible. Because the apparatus does not try to correct the artifacts on the level of 
30 individual projections but processes the whole reconstructed sectional image, it has the 
further advantage to operate quite fast. 
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The sectional image that is processed by the apparatus may be a two-dimensional section 
through an object like the human body, or it may be three-dimensional and represent a 
volume of interest in an object. In both cases, the baseline function is of the same 
dimensionality as the original sectional image. In case of a three-dimensional sectional 
image, the baseline function may be composed of separate two-dimensional baseline 
functions, wherein each of these two-dimensional baseline functions is calculated for a two- 
dimensional slice from the original three-dimensional sectional image. In this case the 
separate two-dimensional baseline functions are preferably fitted to each other in 
neighbouring slices in order to achieve a smooth function all over the whole volume. 

According to a preferred embodiment of the apparatus, the determination of the baseline 
function comprises the following two sub-steps: 

The segmentation of areas from the original sectional image, in which the 
reconstructed X-ray density lies within a given interval. This interval may especially 
comprise the X-ray density of water, such that the segmented areas correspond to 
regions of the object that primarily comprise water. 

Determination of the baseline function based only on data that correspond to said 
segmented areas. Preferably the segmentation is followed by the operation of 
erosion in order to avoid misleading results at the boundaries of the segmented 
areas. 

The restriction of the search for the baseline functions to certain areas of the sectional 
images is especially successful in images of biological objects like the human brain. In this 
case, the image can roughly be subdivided into three classes, i.e. bone, air, and soft tissue, 
the latter corresponding in its properties approximately to those of water. The tissue areas 
may then be taken as a reference region for the desired calculation of the baseline function. 

In a preferred embodiment of the apparatus described above, the baseline function is 
determined by fitting a parametric model function to the data in the segmented areas. The 
model function may be especially a polynomial (preferably of 4th to 6th degree) and/or a 
spline function. 
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According to another embodiment of the apparatus that calculates segmented areas, the 
baseline function is determined by low-pass filtering of data in the segmented areas. This 
approach is especially suited in cases where the segmented areas are spatially distributed 
5 over the whole image area such that the low-pass filtering will be able to fill the gaps in 
between them. 

Still another approach for the determination of the baseline function is based on spectral 
filtering. In this case, the apparatus is adapted to execute the following steps: 

10 

A spectral analysis (e.g. the three-dimensional Fourier decomposition) of the 
sectional images, or, in the case of an apparatus that segments the image, of the 
segmented areas of the sectional image. 

1 5 A composition of the baseline function from only the lower frequency components 

of the resulting spectrum. 

An advantage of such an apparatus is that it may make use of existing algorithms and 
software libraries for spectral processing of images. 

20 

According to another embodiment of the apparatus, image areas outside the object are 
segmented in the original sectional image, and no correction of the original image is done 
there (i.e. the baseline function is set to zero in these areas). By this way stability of the 
results is improved. 

25 

The apparatus may especially comprise a rotational cone beam X-ray device for the 
generation of X-ray projections of an object. These projections may then be used for the 
reconstruction of the sectional image that is further processed by the apparatus. Especially 
in cone beam X-ray projections, there is a considerable contribution of inhomogeneous 
30 artifacts that need to be corrected. 
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The invention further comprises a method for the processing of a sectional image that is 
reconstructed from X-ray projections of an object from different directions, comprising the 
following steps: 

a) The determination of a baseline function that describes spatially slowly varying 
artifacts of the sectional image. 

b) The calculation of a corrected image by compensating the sectional image with the 
help of said baseline function. 

The method comprises in general form the steps that can be executed with an apparatus of 
the kind described above. Therefore, reference is made to the preceding description for 
more information on the details, advantages and improvements of the method. 

In summary, the present invention proposes a retrospective homogenization procedure. In 
a preferred embodiment, the reconstructed volume is segmented into bone, tissue and air 
based on their gray values, and the tissue-regions are used as support in order to fit a 
spatially slowly varying baseline function representing the smooth shape of cupping or 
other inhomogeneities. Finally the inverse of the estimated baseline is subtracted from the 
original slice to correct for the inhomogeneities. 

Beside fast calculation another advantage of the proposed method is that, in contrast to 
projection based correction schemes, the retrospective homogenization does not necessarily 
need detailed knowledge of system or acquisition parameters. 

These and other aspects of the invention will be apparent from and elucidated with 
reference to the embodiment(s) described hereinafter. 

In the following the invention is described by way of example with the help of the 
accompanying drawings in which: 
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Fig. 1 is a diagram of the homogenization algorithm according to the present 



invention; 



Fig. 2-5 depict the processing of a simple example slice in different stages of the 
homogenization algorithm. 

Figure 2 depicts a simple sectional image I through a cylindrical test volume that comprises 
deeply dark and light structures (representing air and bone in a human body, respectively) 
embedded in a "tissue" of a mediate gray value. Moreover, there are some structures of 
small contrast with respect to the "tissue" and a spatially slowly varying inhomogeneity. In 
the lower part of Figure 2, a diagram of the gray value distribution along horizontal line 
through the image (at y = 200) is shown. When a tight gray level window is applied to the 
slice of the reconstructed 3D image in order to find low contrast details, the 
inhomogeneities cause parts of the region of interest to be either below or above the edges 
of the applied window. Low contrast details are therefore much more difficult to find in 
presence of inhomogeneities. 

The basic idea of the present invention is to first find in each slice I of a 3D volume a 
spatially slowly varying 2D baseline representing the smooth shape of cupping or other 
inhomogeneities and second to subtract this baseline from the original image. 
In order to find the baseline, two main problems have to be solved. First, a set of voxels has 
to be determined that can be used as supporting points for the baseline, and second, a good 
fit through these points has to be found. An appropriate selection of supporting points can 
be found by segmenting all voxels into three main classes of matter based on their gray 
value: air, tissue, and bone. Since the human body consists mainly out of tissue with a 
linear attenuation coefficient close to water, this can be done by applying a threshold 
around the linear absorption coefficient of water. The value of the linear absorption 
coefficient of water can be obtained by calculating it based on the applied X-ray spectrum 
or by statistical methods applied to the 3D images. 

To find a fit through all irregularly distributed supporting points, several methods can be 
used e.g. based on low-pass filtering, polynomial, and spline fitting or fitting an arbitrary 
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parametric function. The main requirements are to globally model a smooth shape and to 
not follow local edges. By the choice of the function, application specific prior knowledge 
can be taken into account. Such information could be, e.g., the typical shape, area or gray 
level intensity of the inhomogeneities, that can be derived from previous example fits 
5 specific for the application. Finally, the estimated cupping baseline is subtracted from the 
original slice to correct the inhomogeneities. One of the advantages of the proposed 
retrospective homogenization algorithm is that it does not necessarily need detailed 
knowledge of system or acquisition parameters. 

10 An implementation example of the proposed homogenization algorithm is sketched in 

Figure 1. Images illustrating the different steps of the algorithm are shown in Figs. 2-5 for a 
simple mathematical example slice. 

The algorithm may be carried out by a computer 1 being connected to an imaging 
15 apparatus like a cone beam CT (not shown). The computer will be provided with software 
for executing the algorithm which is illustrated by the flow chart of Figure 1. 

In block 10, the algorithm receives as input the original sectional image I (e.g. a slice of a 
three-dimensional image). By taking into account the cupping of the original reconstructed 

20 volume, a threshold window between -200 HU and +50 HU is applied to the image I in 
block 11. Since at the edges between tissue and bone or tissue and air steep transitions of 
the gray values appear that are pardy within the threshold window but are not beneficial for 
the baseline estimation, these border points are eliminated by a morphological erosion 
algorithm in block 12. Starting from the original slice I, the resulting mask M of voxels 

25 selected as "valid tissue" is thus achieved in block 13 and shown in Figure 3 for the example 
slice of Figure 2. 

A fit through all supporting points of mask M is calculated in block 14 by fitting a two- 
dimensional polynomial of sixth order to the data. The obtained baseline B of block 15 is 
30 shown in Figure 4 for the example slice of Figure 2. Despite the high order of the 
polynomial function, stable results are achieved due to the large number of available 
supporting points. In order to avoid stability problems of the fit outside the imaged object, 
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the region of the object is segmented in block 16 by applying a second threshold followed 
by a four-neighbor connectivity algorithm. Finally, the estimated baseline is subtracted 
from the original slice to correct for the inhomogeneities, and the gray level of the image 
may be shifted to an arbitrary target value like the gray value of water according to the 
5 formula 

I*( X ,y) = I( X> y)-B(x,y) + F 

with p. being a predetermined constant, e.g. representing the gray value of "water". The 
10 result is the corrected image I* in block 17 or Figure 5 for the example slice of Figure 2. 
Note that the homogenization algorithm is not intended to restore exact Hounsfield units. 
Rather, its purpose is to improve the detectability of low contrast details to facilitate the 
application of tight gray level windows. 

15 A similar procedure can also be applied in 3D space. The only difference would be to fit 
the baseline function not for each 2D slice separately, but a joint fit direcdy in 3 
dimensions. Another extension of the described method compares the results of 
neighboring slices demanding for similar functions for the shape of the background. 

20 The estimation of the smooth shape of cupping alternatively can also be done be means of 
a Fourier decomposition and taking into account only the few first Fourier components. 
Depending on the data this can be done without the described segmentation or applied on 
the segmented image M. 

25 Finally it is pointed out that in the present application the term "comprising" does not 
exclude other elements or steps, that "a" or "an" does not exclude a plurality, and that a 
single processor or other unit may fulfill the functions of several means. Moreover, 
reference signs in the claims shall not be construed as limiting their scope. 
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CLAIMS 



1. Apparatus for the processing of a sectional image (I) that is reconstructed from X-ray 
projections of an object from different directions, the apparatus being adapted to execute 
the following steps: 

a) determination of a baseline function (B) that describes spatially slowly varying 
5 artifacts of the sectional image (I); 

b) calculating a corrected image (I*) by compensating the original sectional image (I) 
with the help of said baseline function (B). 

2. Apparatus according to claim 1, characterised in that the sectional image (I) represents a 
10 three-dimensional volume, and that the corresponding three-dimensional baseline 

function (B) is composed of separate two-dimensional baseline functions that are calculated 
for two-dimensional slices of the sectional image (I). 

3. Apparatus according to claim 1, characterised in that the determination of the baseline 
15 function (B) comprises the steps of 

a) segmenting areas (M) from the sectional image (I) in which the reconstructed X-ray 
density lies within a given interval; 

b) determination of the baseline function (B) based only on the data of said segmented 
areas (M), 

20 

4. Apparatus according to claim 3, characterised in that the baseline function (B) is 
determined by fitting a parametric model function to the data in the segmented areas (M), 

5. Apparatus according to claim 4, characterised in that the parametric model function is a 
25 spline function and/or a polynomial, preferably a polynomial of sixth degree. 
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6. Apparatus according to claim 3, characterised in that the baseline function is determined 
by low-pass filtering of the data in the segmented areas (M). 

5 7. Apparatus according to claim 1, characterised in that the baseline function is determined 
by 

a) spectral analysis of the sectional image (I) or the segmented areas (M) of the 
sectional image; 

b) composition of the baseline function from only the lower frequency components of 
1 0 the resulting spectrum. 

8. Apparatus according to claim 1, characterised in that image areas outside the object are 
segmented and excluded from the correction with the baseline function (B). 

15 9. Apparatus according to claim 1, characterised in that it comprises a rotational cone beam 
X-ray device for the generation of X-ray projections of an object. 

10. Method for the processing of a sectional image (I) that is reconstructed from X-ray 
projections of an object from different directions, comprising the following steps: 
20 a) determination of a baseline function (B) that describes spatially slowly varying 
artifacts of the sectional image (I); 
b) calculating a corrected image (I*) by compensating the original sectional image (I) 
with the help of said baseline function (B). 
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ABSTRACT 

Apparatus and method for the processing of sectional images 

The invention relates to an apparatus and a method for the processing of reconstructed 3D 
5 images (I) in C-arm based volume imaging which often exhibit spatially slowly varying 
inhomogeneities caused by inconsistent projection data. To correct the images (I), a 
retrospective homogenization procedure is proposed. The image (I) is segmented (11, 12) 
into principal classes like bone, tissue and air based on their gray values. Only the tissue- 
regions (M) are then used as support in order to fit (14) a spatially slowly varying 2D 
10 baseline (B) representing the smooth shape of cupping or other inhomogeneities. Finally 
the inverse of the estimated 2D baseline is subtracted from the original slice (I) to correct 
for the inhomogeneities. 
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